A new approach to strongly correlated fermion systems: 
the spin-particle-hole coherent-state path integral 
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We describe a new path integral approach to strongly correlated fermion systems, considering 
the Hubbard model as a specific example. Our approach is based on the introduction of spin- 
particle- hole coherent states which generalize the spin-i coherent states by allowing the creation 
of a hole or an additional particle. The action of the fermion system £'[7*, 7; 17] can be expressed 
as a function of two Grassmann variables (7t,7i) describing particles propagating in the lower and 
upper Hubbard bands, and a unit vector field fJ whose dynamics arises from spin fiuctuations. In 
the strong correlation limit, 5[7*, 7; CI] can be truncated to quartic order in the fermionic fields and 
used as the starting point of a strong-coupling diagrammatic expansion in t/U {t being the intersite 
hopping amplitude and U the on-site Coulomb repulsion). We discuss possible applications of this 
formalism and its connection to the t-J model and the spin-fermion model. 



PACS Numbers: 71.10.Fd, 71.27.-fa, 71.30-fh 



I. INTRODUCTION 



Understanding the properties of strongly correlated 
fermion systems remains one of the main goals of 
condensed-matter physics. In narrow-band electron sys- 
tems, the interplay between strong Coulomb repulsion 
(which tends to localize the electrons) and band struc- 
ture effects (which favor their itinerant character) leads 
to a variety of different behaviors ranging from metallic 
Fermi liquid to Mott-Hubbard insulator. Standard weak- 
coupling approaches fail to describe these phenomena so 
that no general theoretical method is available to analyze 
strongly correlated fermion systems. Even for the Hub- 
bard model,Ercl which is supposed to be one of the sim- 
plest (realistic) models of strongly correlated fermions, 
exact solutions or well-controlled approximations exist 
only in a few special cases, lika-in one-dimensionO or in 
the limit of infinite dimensions.El We describe in this pa- 
per a new approach to strongly correlated systems, with 
the (two-dimensional) Hubbard model as primary exam- 
ple. 

Various weak-coupling theories have been applied to 
the Hubbard model and its extensions. The simplest 
of these approaches, the Hartree-Fock theory, supple- 
mented by the random-phase approximation (RPA) for 
the calculation of susceptibilities, correctly predicts an 
antiferromagnetic (AF) ground-state at half-filling in 
both limits of large and small on-site Coulomb repul- 
sion Uu However, it also predicts long-range AF or- 
der at finite temperature in a two-dimensional (2D) sy&j 
tem, in contradiction with Mermin- Wagner theorem.Q 
Many attempts to improve on Hartree-Fock/RPA the- 
ory can be i-fetund in the literature: paramagnon- 
like thepiji£a,Bll3 fluctuation exchange approximaJdaa. 
(FLEX) 00 pseudo-potential Pamuet approach,BEl 
two-particle self-consistent theory.Ej Most of these ap- 
proaches meet with serious difficulties regarding a correct 



description of the physical properties of the 2D Hubbard 
model: absence of long-range order at finite tempera- 
ture, exponential divergence of the magnetic correlation 
length in 1/T at low temperature, existence of Hubbard 
bands and/or precursors of the AF bands in the density 
of states, pseudo-gap at low energy, etc. Besides, irre- 
spective of their success at weak Coulomb repulsion, it 
seems that weak-coupling theories are doomed to fail in 
the strong-coupling regime when the Coulomb repulsion 
exceeds the bandwidth [U ~ ADt with t the intersite hop- 
ping amplitude and D the dimensionality). [See Ref. |lj 
for a detailed discussion of weak-coupling approaches ta 
the Hubbard model.] As emphasized early on by Mott,E£l 
a characteristic property of strongly correlated systems 
is the |existence of local moments already in the metallic 
phaseJj Weak-coupling theories fail to properly describe 
these local moments. For instance, within Hartree-Fock 
theory, local moments are totally absent in the metallic 
phase and appear only in presence of long-range magnetic 
order. 

On the other haadj-the perturbation expansion around 
the atomic limit ,E3~0 which is expected to provide a 
reliable starting point in the strong correlation limit, 
also preseftts its own difhculties in spite of recent 
progress. E^ES It is not completely clear how to handle 
the degeneracy of the ground-state in the atomic limit .l3 
Moreover, the expansion involves two dimensionless pa- 
rameters, t/U and t/T, and therefore breaks down when 
the temperature is much below the electronic bandwidth. 

Nonetheless, it is possible to 4®f¥^ effective Hamil- 
tonian to a given order in t/[/B3~E3 At half- filling, the 
Hubbard model reduces to the Heisenberg model in the 
strong-coupling limit. This model describes local spins 
(i.e. local moments) coupled by short-range exchange 
interactions. Away from half-filling, the effective Hamil- 
tonian in the strong correlation limit is given by the t- 
J model. So far, no satisfying description of a doped 
Mott insulator has emerged from the t-J model. The 
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Schwinger-boson slave-fermion mean-field-theory, which 
provides very good results at halfyfiJlingJ^El has not been 
as successful in the doped case.tM^ On the other hand, 
the slave-boson mean-field theorjH does not even repro- 
duce the known results at half-filling, so that its predic- 
tions (near half-filling) are highly questionable. 

A very appealing path integral formulation of the t-J 
model, based on the spin- hole coherent-states. jJl, C), has 
been introduced by Auerbach and Larson.Efl The exis- 
tence of local moments is obvious in this formulation, 
since a singly occupied site r is described by a spin- 
coherent state I Or). The set of Grassmann variables 
C = {Cr} allows for the presence of holes. The spin- 
hole coherent-state path integral was used to study spin 
polarons in the (large-S*) semiclassical limit of the t-J 
model. 

Schulz has derived a path integral formulation of the 
Hubbard modeL-wthich also exemplifies the existence of 
local moments.L3'Ej The Coulomb repulsion is treated 
within a large-C/ Hartree-Fock approximation, whereas 
the SU(2) spin-rotation symmetry is maintained by intro- 
ducing a fluctuating spin-quantization axis in the func- 
tional integral. The resulting effective action S'[7*,7;ll] 
is expressed in terms of two Grassmann variables (7|,7x) 
which describe particles propagating in the lower (LHB) 
aniijipper (UHB) Hubbard bands, and a unit vector field 
O.EZI A singly occupied site r corresponds to a local mo- 
ment pointing in the direction of Similar ideas, in 
view of Monte Carlo simulation, iave been recently dis- 
cussed by Bickers and Scalapino.E3 

In a previous publication, Pairault and the present 
author have reporte(L-|a systematic t/U expansion of 
the Hubbard model.t3 As in Schulz's work, central 
to this approach is the introduction of a fluctuating 
spin-quantization axis. By adapting the method of 
Refs. the t/U expansion was generated by a 

Grassmannian Hubbard-Stratonovich transformation of 
the hopping term. As in Ref. the effective action 
5[7*,7;r2] of the Hubbard model in the strong-coupling 
regime is expressed as a function of two Grassmann vari- 
ables and a unit vector field. We found however that the 
Hartree-Fock treatment of the Coulomb repulsion is not 
sufficient in the strong correlation limit, since it misses 
processes of order t/U . In particular, it does not allow ta 
recover the t-J model when the UHB is integrated out.t£l 

The aim of this paper is to further develop the strong- 
coupling expansion introduced in Ref. |3^. (i) The func- 
tional integral formulation is derived from a completely 
different perspective, starting from spin-particle- hole co- 
herent states (Sec. ||). The latter generalize the spin-i 
coherent states by allowing the creation of a hole or the 
introduction of an additional particle. This new deriva- 
tion emphasizes the connection with the spin coherent- 
state path integral formulation of the Heisenberg model, 
(ii) As in Ref. the strong-coupling expansion is gen- 
erated by two successive Hubbard-Stratonovich trans- 
formations of the intersite hopping term. The result- 
ing action S'[7*,7;r2] is obtained to all orders in t/U 



(Sec. [HA). This action contains interaction terms to 
all orders which are determined by the (exact) atomic 
vertices. In the strong correlation limit, 5[7*,7;J1] can 
be tr uncat ed to quartic order in the fermionic fields 
(Sec. HIC). The action is then entirely determined 
by the single-particle atomic Green's function and the 
two-particle atomic vertex, (iii) We discuss in detail 
the strong-coupling perturbative expansion and briefiy 
point out possible applications of our formalism ■wjhjfhi 
are develo ped in detail in separate publication^^ 
(Sec. pID| ). Finally, we discuss the connection of this 
formalism to the t-J model and the spin-fcrmion model. 



II. SPIN-PARTICLE-HOLE COHERENT-STATE 
PATH INTEGRAL 

In this section, we define the spin-particle-hole coher- 
ent states and derive a path integral representation of 
the partition function. To be specific, we consider the 
Hubbard model defined by 



H^-t ^ {cl„Cr'a+'i^-C.) + U\ 
(r,r') ,a 



[flr^flri, (2.1) 



where Cro- is a fermionic operator for a cr-spin particle 
at site r (cr =t,i), fira = cl^Cra, and (r, r') denotes 
nearest neighbors. For simplicity, we consider a bipartite 
ZJ-dimensional lattice. We denote by fj, the chemical po- 
tential and consider only hole doping (i.e. /i < U/2). We 
take h — ks — ^ throughout the paper. 



A. Spin-particle-hole coherent-states 

1. Definition 

We first consider a single site. A basis of the Hilbert 
space is 6 = {|0), | T>, | i>, | Ti)}, where |0), | T) = c\\0), 
I J,) = c||0) and | ti) — are the empty, singly 

occupied (with spin up and spin down) and doubly occu- 
pied states, respectively. Instead of using the states | t) 
and I I) to describe a singly occupied site, we introduce 
the spin-i coherent state 

= cos^e-5(v+V')| |) + sin^e5(^-^)| |), (2.2) 

where 9 and are the polar angles determining the direc- 
tion of the unit vector Vi(0, ip). Thevchoice of ip is free and 
corresponds to a "gauge" freedom.CJ B'^ = {|0), jJl), | ti 
) } will be used to construct the spin-particle- hole coher- 
ent states. Although it is not (for fixed ft) a basis of the 
Hilbert space, we will show how the whole Hilbert space 
can be described by varying ft. 
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We now enlarge the physical Hilbert space by intro- 
ducing three bosonic operators (e, p and d) and rewrite 
the states of B'^ as 



|0) = et|vac) =7t|JI), 
|f2)=pt/t|vac), 

n>=dV7/iVac) = -7l|f^), 



(2.3) 



where |vac) denotes the vacuum of the enlarged Hilbert 
space. The latter contains unphysical states which can 
be eliminated by imposing the constraints 



q\ '^ ^ f\ /t 



: 0. 



1=0, 
= 0, 



(2.4) 



The (slave) bosonic particles are similar to those intro- 
duced by Kotliar and Ruckensteira although here the p 
boson is spinless. Once we have introduced the spin- 
i coherent state we need two operators, 7| and 

7l, to allow for the creation of a hole or an additional 
particle. These two operators are Hubbard operators^ 
(7-f = |0)(f2| and7x = — |r2)(t| |), and have the following 
expression, 



7T = e^t. 



7i 



(2.5) 



in terms of the slave bosons and the fermionic operators 
/ct. Note that all the operators are defined with respect 
to the spin-quantization axis Jl, so that a /|- (or 7|-) 
particle has its spin pointing along H. Thus, the LHB 
is populated only with up-spin particles (in the spin ref- 
erence frame defined by fJ). This also implies that only 
down-spin particles can be introduced in the UHB. The 
operators 70- play a fundamental role in our approach. 
Although they app ear at this stage as composite oper- 
ators [see Eqs. (2.5)], we shall show in Sec. [H how the 
partition function can be written as a functional integral 
expressed in terms of the elementary fields 70- and the 
unit vector field fJ. 

The spin-particle-hole coherent states are defined by 

\n, C) = exp(eet + ppt + rfjt _ /^/t _ fJl)\^^c), (2.6) 

where f-\,fi are Grassmann variables and e,p,d c- 
numbers. We use the short-hand notation C = 
(f^l , fl,e,p,d). By using the constraints (2.4), we can 
rewrite the spin-particle-hole coherent states as 



\nx)^e\o)-pf^\n)-df^f^\n)- 



(2.7) 



Eqs. (p.q) and (2.7) extend the definition of the spin- 



coherent states by allowing for the presence of a hole or 
an additional particle. 



2. Scalar product 

We consider the scalar product between the states 
|r2,C) and \rt',('). If = i7' all fermionic variables 
are defined with respect to the same quantization axis. 
We then have the standard resultcS 

= exp(e*e' +pV + d*d' + ^ /;/;). (2.8) 



When O 7^ Jl', we use Eq. (2.7) to obtain 

{n, C|0', C') = e*e' +p*f^f\p' + d*d'flf*f\f[ 

+p*f^f\p'i{n\n') - 1). 



(2.9) 



For rj = r2', this expression should be equivalent to 
Eq. (^, so that 

exp(C*C') = e*e'+p*flf{p' + d*d'flf\f\f[ (2.10) 



provided the constraints ( |2.4D are satisfied. Similarly, by 
keeping only the terms which are compatible with the 
constraints, we obtain 

exp(CC' + f;p*f\p'{W) - 1)) - e*e' 

+p*f^f\p'{n\n') + d*d'fir^f\fl (2.11) 

a result vali d wh ateve r the value of Vt and VI' . Com- 
paring Eqs. (|2j) and ( ^.ll[ ), we arrive at the following 
expression for the scalar product of two spin-particle- hole 
coherent states in the physical Hilbert space: 

{n,c\^'X') = exp(cc' + /TV/|p'((f^|f^') - 1)). 

(2.12) 

3. Resolution of the identity 
We seek a resolution of the identity in the form 

^1^1 dCdCe-"^'^^'\n,0{nx\P^i, (2.13) 

where / is the unit operator and 

a|CP = aelep + \p\' + + ^ (2.14) 

cr 

J\f, Q!e a nd a d are constants to be determined in order 
for Eq. ( |2.13| ) to be satisfied. The projection operator 



(2.15) 



ensures that the resolution of the identity acts only in 
the physical Hilbert space. The measure in Eq. (2.13) is 
defined by 
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dcdc = '4^'4^'4^df^df,df:dn, (2.16) 



where 



2iTT 2in 2i7r 



dz*dz 1 

= —dRe(z)dlm(z) 



(2.17) 



for a complex variable z. Integrating over ^, we rewrite 
Eq. {WJ^ as 



i^^Xf^i , in)(n 



= i. (2.18) 



Since the spin coherent states satisfy 2{An)^^ J dfl\fl) (r2| : 
I T)(T I + I I (resolution of the identity in th e sp ace 
of singly occupied states) , we conclude that Eq. ( p.l3| ) is 
satisfied for 



M = 8, 

ae = oid = 2. 



(2.19) 



The constants ae = ad (7^ 1) are necessary to avoid an 
over-counting of the empty and doubly occupied states: 
for a given fJ, B'f^ already contains the states |0) and 
I ti)- Varying Q, with ae = ad = 1 would lead to an 
over-counting of these states. 

Note that Eq. ( 2.13| ) bears some similarities with the 
resolution of the identity in the Hilbert space with no 
double-0|Geupancy obtained from the spin-hole coherent 



states.l 



4-. Trace 



The resolution of the identity ( p. 13 ) allows to com- 
pute the trace of a given operator O in terms of the 
spin-particle-hole coherent states. The derivation is 
standardly and one obtains 



B. Path integral for a single site 

In order to express the partition function as a path in- 
tegral, we use Eq. (2.20) and divide the "time" interval 
/3 into M steps: 



47r 



X (n, ^\Pe-<"~^''^^ ■ ■ ■ e"'^"-^'^\n, Q, (2.22) 



where e = j3/M. We the n intr oduce (M — 1) times the 
resolution of the identity ( 2.13|). The details of this pro- 
cedure are given in Appendix lAl. One finds 

n '^^dQdCudX^ 




fc=i 



47r 



M 

xexp<j^[-aiap + aa-i 



+/ffePfc/Tfe-iPfe-i((^fel^fc-i) ^ 1) ^ <^Kk,k-i 



where 



k,k-l — 



{nkXk\Kk\nk-i,Qk-i) 



(2.23) 



(2.24) 



kk = H - iiN + iX'-k^ Q^k^ 



Afe = (X^k\x!'^k) denotes Lagrange multipliers which im- 
pose the constraints (2.4). 

An important difference with the standard path inte- 
gral for a system of fermions and bosons comes from the 
presence of the constants ae = a^ (7^ 1) which ensure 
that the empty and doubly occupied states are not over- 
counted. Taking the continuum time limit (e — > 0) from 
Eq. ( p. 23 ) would then lead to an infinite chemical poten- 
tial ^e — —{oie — l)/e for the e boson: 



^^•^^^/l^ / dCdCe-^^^\nX\Pd\n,0, (2.20) Y.el{ek-.-a.ek)^-[dre*{^ 



where C = (-/T,-/i,e,p,d). 

As a simple application of Eq. ( 2.20| ), we calculate the 
partition function for a single site. Imposing the con- 
straints (2^) by using Eq. (2/?), we easily obtain 



Ot H e, 



(2.26) 



-(ad 



Zat = Tre 



-p(H-tj,N) 



1 



1 



„/3(2p-C/) 



^e"d 



(2.21) 



where N is the total number of particles and (3 — \/T is 
the inverse temperature. 



and a similar result for the d boson (i.e. ^.d 

i)A)- 

At half- filling, fj, = U/2 due to particle- hole symme- 
try. The contribution of the empty and doubly occupied 
states to the partition function is exponentially small at 
low temperature (T ^ U/2) and can be neglected. It 
is then possible to take ae = ad = 1, and the contin- 
uum time limit is well defined. This argument fails away 
from half-filling in the single-site case. However, we are 
interested in the full Hubbard model, where the chemi- 
cal potential is determined by the filling of the Hubbard 
bands in the strong correlation limit. In presence of a 
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small concentration of holes, fi ~ W/2 is located near 
the top of the LHB (of total width W). In the low- 
temperature limit T ^ W, the system in the atomic 
limit is in its ground-state with exactly one particle per 
site. The contribution of the empty and doubly occupied 
atomic states to the partition function is again negligible, 
which allows to take ae — ad — 1- This argument is par- 
ticularly clear in the strong-coupling expansion (Sec. |l|) 
where the atomic action is used only in the calculation of 
the connected atomic functions G^'^. The latter are ob- 
tained from the atomic limit {t = 0) of the full action, the 
chemical potential being determined by the full Hubbard 
model. 

Taking ag = = i-iand ignoring the overall normal- 
ization constant A/'*^Ji3 we write the partition function 
as 



exp-! 



E 

fc=i 



[CkiCk-Ck-i) 



-f*kP*khk-iPk-ii{^k\^k^i) - 1) + eKk^k-1 



(2.27) 

The continuum time limit can now be taken without dif- 
ficulty (see Appendix^: 



Zat ^ JvflJvxJ V[f,e,p,d]( 
where the action is given by 



(2.28) 



c _ o(0) , 
"^at — Oat + 



dTf*p*f^p{n\n) 



c(0) _ 
"-"at — 



dr 



a 

d*{dr + U + tX^^'^ -iY,>^a^)d 



+e*{dr + iX^^'^)e+p*{dr + ?;A(i) - iXf^)p 



(2.29) 



S'at^ is similar to the actian obtained from the Kotliar- 
Ruckenstein slave bosonaij with p-^ = p and without the 

Pl boson. The difference between S'at and sl^\ which 
arises frorn-tJic dynamics of $7, is nothing but the Berry 
phase termlia = (r2|r2) of a spin- i = dr\n)). The 
Berry phase term is alive whenever the state is singly oc- 
cupied (/*/| = 1 and p*p =1). 



C. Path integral for the Hubbard model 



where the unit vector is now site dependent [Tt = 
{fir} and C = {Cr}]- The fermionic operator f^a is de- 
fined with respect to the local spin-quantization axis fir- 

The procedure to write the partition function as a path 
integral is similar to the single-site case. The only con- 
tribution which is not purely atomic comes from the in- 
tersite hopping term. The latter modif ies K k,k-i (or, in 
the continuum time limit, Kk^k) [Eq- (2.24)]. As shown 
in Appendix |^ 

{Ctk, Ck\clcr'\flk, Ck) = llkRlk^r'kJr'ki^k, Ck\^kXk), 

(2.31) 



where Cr — (CrtjCrx)"^- 



(2.32) 



is a SU(2)/U(1) matrix which rotates the spin- 
quantization axis from z to n{9rk,y^rk)- It satisfies 
RrkC^zRlk ~ ^rk ' c , whcrc CT — {(Jx, (7y, fJz) stands for 
the Pauli matrices, jrk — (7rTfc; TrXfe)'^, and 



7rTfc = KkPrkfr-fk, 
7rifc = Vlkdrkfrlk- 

The action of the Hubbard model then reads 



(2.33) 



7r', 



(2.34) 



where S'at is the genera lizati on to N sites of the atomic 
action given by Eqs. ( ^.29 ). In the Hubbard model, 
trr' equals t if r and r' are nearest n eigh bors and van- 
ishes otherwise. Since the constraints (2.4) are preserved 
under the time evolution determined by the action S 
[Eq. (2.34)], we can replace the functional integral over 
A(t) by an integral over a set A = (Ar^^Ar^) of time- 
independent Lagrange multipliers: 



VX 



.A.ni 



2-K 



n 



2-K 



(2.35) 



The definition of the spin-particle-hole coherent states 
can be extended to the lattice case: 



J7, = exp(^^(ere]: + Prp\ + drd\ - ^ 



I vac), 



(2.30) 



Equations (2.29) and (2.34) provide a spin-rotation- 
invariajHt slave-boson formulation of the Hubbard 
model.Ej They have been derived recently without refer- 
ing explicitly to the spin-particle-hole coherent states, 
but by introducing a fluctiiating spin-quantization axis 
in the functional integralHj 
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III. STRONG-COUPLING EXPANSION 



where 



In this section, we recast the spin-particle-hole 
coherent-state path integral in a form suitable for a per- 
turbative expansion with respect iot/U . This is achieved 
by performing two successi ve tra nsfo rmati ons of the in- 
tersite hopping term (Sees. [II A and p^II B ). The result- 
ing action is expressed as a function of two fermionic 
fields (71 and 74 ) which describe the propagation of par- 
ticles i n the LHB and UHB, and a unit vector field rj. 
In Sec. |IIIC , we show that this action can be truncated 
to quartic order in the fermionic fields in the strong cor- 
relation limit. Finally, in Sec. HID we show how this 
effective strong-coupling action can be used as the start- 
ing point of a diagrammatic perturbative expansion with 
respect to t/U , and discuss possible applications of our 
formalism. 



Z^,[U\ = JdxJ V[f,e,p,d]t 



Z. 



(3.5) 



is the "partition function" in the atomic limit {t = 0) for 
a given configuration of n. zj^^ is the partition func- 
tion obtained from sl^\ We emphasize here that the 
integration of the Lagrange multipliers and the bosonic 
fields has been done exactly, ^^[n] is the action of the 
spin degrees of freedom in the atomic limit. Treating 
A'^ = (r2r|i^r) in perturbation, we obtain to lowest order 
(in a cumulant expansion)Ej 



A. Grassmannian Hubbard-Stratonovich 
transformations 

Following Refs.@0@ we decouple the intersite hop- 
ping term in Eq. ( |2.34 ) by means of a Grassmannian 
Hubbard-Stratonovich transformation. The partition 
function becomes 

Z = J -Dndetii) J I?Mcxp{-^^:f-,Vb} 

a, 6 

X dX I?[/,e,p,d]exp{-5at + 5Z(V'a7a + c.c.)}, 

(3.1) 

where if) is an auxiliary fermionic field whi ch couples to 
the (composite) field 7 defined in Eq. ( 2.32 ) . We use the 
notation 



r 



(3.6) 



Sb is a collection of Berry phas e te rms for spins local- 
ized at the lattice sites. In Eq. (iA), W[ip*,'il};fl\ is the 
generating functional of the connected atomic Green's 
functions 



(3.7) 



Note that the G^'^'s are defined for a given configuration 
of the field n. 

The partition function is thus written as 



Z = Z, 



(0) 



Vndet{i) V[iP]e 



and denote by 



^ „ J a.b 



(3.3) 



the rj-dependent intersite hopping matrix. The 
SU(2)/U(1) matrix rotates the spin-quantization axis 
from z to Jlr [see Eq. ( ^.32| )]. The diagonal elements of 
trr' correspond to intraband propagation, while its off- 
diagonal elements describe transitions between the LHB 
and UHB. Note that det(t) is a function of Vl and should 
therefore be kept explicitly in the functional integral. 

Performing the functional integral over the fields 
/, e,p, d and the Lagrange multipliers A, we obtain 

I d\( I?[/,e,p,d]exp{-5at+5](^:7a+c.c.)} 

= Zat[n]e'^['^'''^'"l (3.4) 



"0; rj] can be obtained explicitly by inverting 
Eq. (Q: 



(i?!)2 

R=l ^ ' ai,bi 



(3.9) 



The primed summation in (|3.9D reminds us that all the 
fields in a given product tp^i ' ' ' "061 share the same value 
of the site index. 

In order to completely determine the action 
S[ip* we need to compute the atomic Green's 

functions G^'^ from the action S^f The procedure to 
obtain these quantities can be found in Appendix ^ and 
in Ref. Here we only quote the final results. The 
single-particle Green's function is 



6 



= Gi")-i-sgn(a)AO, (3.10) 
l(--l) for cr =t (i)- Gf") is the atomic 



where sgn(o-) 

Green's function corresponding to S'^^j . 
perature hmit 



:.(o) 



In the fow tern- 



where rj 

;(0), 



GV'ituj) 



(r)=0(-r + 7y)e^^ 
(r) = -9{t ~ T^)e^t'-^^\ 
0+. In Fourier space, Eqs. 



(3.11) 



and G^^Hiuj) 



'I 



( 3.11 ) become 



(lU! + fj, 

where uj is a fe rmionic Matsubara frequency. As shown 
in Sec. (IIIC), the i?-particle Green's functions G^'^ 
[R > 2) arc involved only in (virtual) interband tran- 
sitions. The typical energy scale (~ U) for the latter is 
much larger than the typical energy scale for spin fluc- 
tuations J — At'^/U) so that the effect of the Berry 
phase term A" can be neglected (i.e. we can use S^^^ 
instead of S^t to compute G^'^ for R> 2). We find 



G 



Rc 



(i? > 2), 



(3.12) 



which implies that there is no intraband interaction for 
the ip field. The two-particle atomic Green's function 
G^^'^ is thus entirely determined by G^j_| , or equivalently 



by the two-particle vertex (see Eq. (C12)) 



0^3 is fixed by energy conservation. 
IlII C that the effective action in the 



r||_ll(icji, 1^2; jti;3(ia;4)) ~ iuj2 — — U, (3.13) 

where UJ4 = uJi + uj2 - 

We shall show in Sec. 

strong-coupling regime {t/U <^ 1) is entirely determined 
by the knowledge of G and G^^"^ (or F^^) so that we do 
not require the knowledge of G^'^ for i? > 3. 

We could now follow Ref. |2^ and use the action 
S[il)* to perform a diagrammatic perturbative ex- 
pansion with respect to the intersite hopping amplitude 
t. The matrix structure of trr' distinguishes processes 
that do conserve the number of doubly occupied sites 
from those that involve interband transitions, thus al- 
lowing for a t/U expansion. [In Ref. the auxiliary 
field tp^ couples to the original fermionic field c^. The 
resulting perturbative expansion involves two dimension- 
less parameters, t/U and t/T, and breaks down at low 
temperature.] However, the drawback of using the ac- 
tion S[ip*,tlj; Jl] is that the field ip has no direct physical 
meaning. It is rather the fields to which ip-^ and cou- 
ple, namely the composite fields e*pf-\ and p*dfi, which 
have a direct physical interpretation in terms of parti- 
cles propagating in the LHB and UHB. This suggests to 
perform a second Hubbard- Stratonovich transformation 
of the intersite hopping term. Denoting by 7 the auxil- 
iary field of this transformation, we rewrite the partition 
function as 

Z = Zj°^ JvnJ 2?[7, V] exp{-5B[f^] + ^*Jablb 



+ I](7a>a + C.C.) + W[r,i^; ^]}- 



(3.14) 



Integrating out the i/j field, we obtain 



/ V[j]e^p[Y,i:iabJb + W[j*,r,n]}, (3.15) 



Note that 7 is now an elementary field. [We expect 
7 to have the same physical meaning as the compos- 
ite fields defined in Eq. (2.33), hence the same nota- 
tion.] VF[7*, 7; fi] is the generating functional of the con- 
nected Green's functions G^'^ calculated with the action 
— W[^*,ip; ft]. Since the latter is local (in space), so are 
the G^^'s. Z[n] is defined by 



Z[n] = J I?[V;]e^[^*''^'"] 
-det(-G)e-^'["l. 



(3.16) 



det(— G) comes from the Gaussian part, J2a b^at'tpl'tpb, 
of the action —W[il!*,ip;n]. Non-Gaussian terms give 
the contribution S"[r2]: 



oc 



=-s'[n] _ 



= ( exp' 



V (-1) 

^ (i?!)2 ^ 



(3.17) 



where the average is to be taken with the action 
J2a bGab4'l''Pb- 5" [J7] Can be represented as a sum of 
Feynnian diagrams, the bare propagator being —G^^ 
and the vertices the atomic Green's functions G^'^ (with 
R>2). An example of diagram is given in Fig. ^ [The 
symbols used in the Feynman diagrams are defined in 
Fig. |l|.] This diagram, as well as the other diagrams con- 
tributing to S"[ll] does not have any physical meaning. 
As will be shown below, S' [Jl] does not contribute to the 
final result. "Physical" and "anomalous" (i.e. without 
"physical" meaning) diagrams can be defined rigorously 
by considering the per turb ation theory based on the ac- 
tion Sitp*,^};^.] [Eq. dU)]. Within the latter, all dia- 
grams can be represented by atomic Green's functions 
connected by the intersite hopping matrix t. The per- 
turbation theory based on the action S'[7*,7;r2] gener- 
ates diagrams that cannot be represented in this way. 
These "anomalous" diagrams should not contribute to 
physical quantities. Indeed, we will find that the total 
cont ributi on of anomalous diagrams always vanishes (see 
Sec.|IHBD. 
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G- 



G„ 




G- 




G-1 

FIG. 3. A Feynman diagram contributing to F. 



t = 



FIG. 1. Definition of the various symbols appearing in 
the Feynman diagrams. Connected atomic Green's functions 
are represented by a filled circle and atomic vertices by poly- 
gons, with the exception of which is also shown as a 
filled circle. The intersite hopping matrix t is represent ed by 
a dashed line, and the propagator Q^^"^ defined in Eq. (3.37) 
by an empty circle. 



A result similar to Eq. ( 3.19| ) holds for higher- 
order Green's functions. The simplest diagram 
for the two-particle Green's function Gl^'l u h — 
(V'aiV'asV'baV'bi)-^ (shown in Fig. |a) gives the contri- 
bution -T^y^ t^b2' where r]^,a2Mb. two-particle 



vertex. All other diagrams are anomalous (Fig. In 
the same way, G^^^'^ can be written as the sum of F^^^ and 
an anomalous contribution f (Fig. ||). More generally, 
we find 



Q- 



FIG. 2. A Feynman diagram contributing to 5'[f2]. 
Noting that det(-G) = we obtain 

s[r,r,^] = s'[n] ^Y.^%ab - wh*,i;^]- (3.i8) 



6.} 



f J?. 



ft.} 



{R > 2), (3.20) 



In order to completely determine the action S['y*,^;d], 
we need to compute the (local) Gr een 's functions G^'^ 
from the action -^^[1/^*, -0; fJ] [Eq. 1^^]. 

Let us first consider the single-particle Green's func- 
tion Gab = —{^aipD-w- Retaining only the Gaussian 
part of the action —W, i.e. J2a bGab'ipl^b, one obtains 
G = — G^^. More generally, we have 



Gab - -G^^ 



Tab 



(3.19) 



where F^ is the atomic i?-particle vertex and f ^ denotes 
anomalous contributions. 



(a) 



G- 



G- 



X 



G- 



where f can be represented as a sum of Feynman dia- 
grams arising from the vertices G^"^ (with R > 2). An 
example is shown in Fig. ^. It is easy to see that all 
the diagrams contributing to F are anomalous. The fact 
that F has no physical meaning can be understood as fol- 
lows . We can view F as a self-energy correction to G [see 
Eq. ( 3.19| )]. Since G is the exact atomic Green's function, 
we do not expect any local (i.e. atomic) correction. We 
will see in Sec. Ill B| that the role of F is to cancel other 



anomalous contributions arising in the strong-coupling 
perturbative expansion. 




(b) 



FIG. 4. Feynman diagrams for G^^". (a) -F". (b) Anoma- 
lous diagram. 
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(a) 




(b) 



FIG. 5. Examples of Feynman diagrams contributing to 
(a) Contributions to F^''. (b) A contribution to the 



anomalous part F 
legs denote . 



All the filled circles with two external 



We can then rewrite the action as 



(3.21) 

Note that this result has been obtained without any 
approximation and gives the exact value of the parti- 
tion function if one carries out the functional integral 
over the fields 7o- and SI. By performing two successive 
Hubbard- Stratonovich transformations, we have summed 
the atomic contributions to the single-particle propaga- 
tor and the i?-particle vertices (i? > 2). As a result, 
the action S[j*,j;fi\ is essentially parametrized by the 
intersite hopping matrix t and th e ato mic vertices G"^ 
and (i? > 2). [We show in Sec. |IIIB| how to deal with 
the anomalous vertices T^.] Using S'[7*,7;r2] instead 



of 5'[/,e,p,d,A;fl] [Eq. (|2.34D ] or S[ip* ^ijj^n] [Eq. 
presents two main advantages. First we now deal with 
only three fields whose physical meaning is clear. 7| 
and 7x describe particles propagating in the LHB and 
UHB, respectively, while the dynamics of fl arises from 
spin fluctuations. [Note however that the Green's func- 
tion of the original fermions (those involved in the deflni- 
tion of the Hubbar d mo del) is not merely the propagator 
— (7a 7b) (see Sec. [IHD ).] Second, 5[7*,7;r2] turns out 
to be a convenient starting point in the strong correla- 
tion limit t/U ^ 1. Indeed, for ?7 ^ oo (no interband 
transition), only the Gaussian part of the action needs 
to be considered. The lowest-order corrections in t/U 
are determined by a small number of atomic vertices. In 
practice, S'[7*,7; $7] can be truncated to quartic order in 
the fermionic fields and is then parametrized only by 
and F" (Sec. pTcl). 



B. Anomalous diagrams 

In order to understand the role of anomalous diagrams, 
it is useful to fir st consider the atomic limit {t = 0) of 
the action ( |3.2l| ). In that limit, one has to recover the 
known (exact) results. 

Let us consider first the partition function. Integrating 
out the 7 field, one obtains 

Zat[f2]=det(-G-i)e-^'["l-^"["l, 



e-^"["l = (exp 



a,b 



oo 



■ • • 76i 



(3.22) 



where the average is to be taken with the Gaussian action 
~ b^aGabJb- In Order to obtain the correct result 

= Zl°^y"l?fle-^«["1, (3.23) 

one must have S'[rt] + S"'[J7] = 0. This result can be 
obtained by inspection. Consider for instance the dia- 
grams containing once and only once the vertex F^^ (and 
no higher-order vertex). The diagram contributing to 
S'[fl\ is shown in Fig. g and has the overall factor —1/2. 
There are two contributions to S"'[r2]. They come from 
the terms (j^^ablb) and (Taia2Mb2^a^^a2^b2lb,) where 
the average is taken with the action — ^ lt,G~^^b [see 
Eq. ( 3.22 )]. Both contributions correspond to the dia- 
gram of Fig. H but with the factors 1 and —1/2. We 
conclude that the total contribution to S' + S" vanishes. 

In the same way, we can calculate the single-particle 
Green's function —{'^all) (^°^ ^ given configuration of 
n) and verify that we do obtain the correct res ult in the 
atomic limit, namely Gab- From the action (B.21), we 
read off the inverse (atomic) propagator G^^ — F. Self- 
energy corrections due to vertices i? > 2 should therefore 
cancel F. Again this result can be proved by inspec- 
tion. Consider again the diagrams containing once the 
vertex F^^. The contribution to F is shown in Fig. |^. 
The one-l oop self-energy correction coming from F^^ (see 
Eq. ( 3.21 )) produces the same diagram but with opposite 
sign. 

Thus we come to the conclusion that anomalous dia- 
grams cancel in the atomic limit. We expect this prop- 
erty to hold also when t ^ Q, since anomalous diagrams 
do not have any physical meaning. Consider for instance 
the anomalous diagrams shown in Fig. ^ which con- 
tribute to the effective action of the spin degrees of free- 
dom S'[J7] = — lnZ[r2]. [Fig. ||b shows the correspond- 
ing "physical" diagram.] Working out signs and symme- 
try factors, we find that the sum of the two diagrams of 
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Fig. ^ vanishes if f is approximated by the contribution 
shown diagrammatically in Fig. ^. 



(a) 





must necessary contain closed loops of 7^ particle. In 
the absence of particles in the UHB, these loops vanish 
(see Appendix ^). As a result, the effective action of the 
LHB is Gaussian to leading order in t/U: 



r •' 

dTlr-\ir],v']lr']- (3.26) 



(b) 



FIG. 6. (a) Two anomalous diagrams for the free energy. 
Their sum vanishes when V is approximated by the contribu- 
tion shown diagrammatically in Fig. |^. (b) The corresponding 
"physical" diagram. 

As expected on physical grounds, anomalous diagrams 
always cancel in the calculation of physical quantities. 
We can then work with the action 



E 

a.b 

CO 



ETpiwErfa^M^-'-T.. (3-24) 



(i?!)2 

R=2 ^ ' ai,b, 

provided that we discard all anomalous diagrams. 



C. Effective action of the Hubbard model in the 
strong correlation limit 

In this section, we show that only the quadratic and 
quartic parts of the action S'[7*,7;J1] need to be con- 
sidered in the strong-coupling limit. This is achieved by 
considering the effective action S'lhb [7* , 7t ! of par- 
ticles propagating in the LHB. To hrst order in t/U, 
S'lhb [7|, 7t; ^] should correspond to the t-J model. This 
result will guide us in the determination of the leading 
terms in the action of the Hubbard model in the strong 
correlation limit. We emphasize however that our pri- 
mary goal is not to determine the effective action of car- 
riers in the LHB, but to derive an effective action for the 
full Hubbard model in the strong correlation limit. As 
discussed at the end of this section, there are advantages 
in working with the Hubbard model instead of consider- 
ing the t-J model obtained by integrating out the WIB. 

S'lhb is obtained by integrating out the field:E3 



3[7|,7t;^^] = 



(3.25) 



glcctcd. Since F^'^^..^. = 0, diagrams contributing to Sj 



(0) 
LHB 



The action (3.26) describes particles propagating in the 



LHB and interacting with spin fluctuations via the gauge 
field A'^ and the rj-dependent intersite hopping matrix t 

(i.e. Sl^b describes fermions coupled to a U(l) gauge 
field). A detailed discussion of its physical meaning can 
be found in Sec. III.B.l of Ref. ||. 

We now cons ider the contributions of order t/U to 
S'lhb [Eq- ( |3.25| )] . They are represented by the diagrams 
of Fig. 0a. Each of these diagrams contains two inter- 
band transitions. The vanishing of closed loops with 
particles ensures that all the diagrams to be considered 
are of the type shown in Fig. |^a. Examples of vanishing 
diagrams are shown in Fig. 



G 



i t 



(a) 





(b) 



At zeroth order in t/U, interband transitions are ne- duce7-f(r') 



FIG. 7. (a) Contributions of order t/U to the effective ac- 
tion S'lhb. (b) An example of vanishing diagram. Slashed 
dashed lines indicate interband transitions. 

Let us first focus on the quadratic and quartic contri- 
butions generated by the integration of the UHB (first 
two diagrams of Fig. ^). These terms are non-local in 
time. However, since they involve (virtual) interband 
transitions, they can be approximated by local vertices. 
From the equations of motiortiii the atomic limit, we de- 



-11(7 



'7T 



(r) E3 Using the latter equation 



and integrating over the time differences, the non-local 
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vertices are then approximated by local vertices. For in- 
stance, the two-point vertex (first diagram of Fig. Qa) is 
approximated as 

J2 J drdr'7;T(r)irT,r'i(r)Gi"Hr-r')4u,r"T(^')7r"T(^') 

r,r' ,r" 

--JjYl J drj*^{r)irhr'l{T)tr'l,r"^{Thr"-f{T), (3.27) 



where we have used Eq. ( 3.11 ). Since we can ignore the 
effect of spin fluctuations on (virtual) interband tran- 
sitions, we have replaced by cj^''. If we also ap- 
proximate the quartic term by a local vertex and ne- 
glect higher-order vertices, we end up with an effective 
action which is nothing but the action of the t-J model 
(including the so-called pair-hopping terms) in the spin- 
hole coherent-state path integralH [See Ref. ^ for a de- 
tailed derivation.] It is crucial here that the effective 
action is local in time, otherwise it would not derive 
from an Hamiltonian. It should be noted that a simi- 
lar approximation is made in the standard derivation of 
the t-J model, which starts from the resolvant operator 
R{E) = [E—H)^^. Projecting out states with double oc- 
cupancy yields an effective "Hamiltonian" Heg{E) which 
depends on the energy E. The Hamiltonian of the t-J 
model is obtained by retaining terms of order 0(t/U) and 



replacing Ei^W the energy Eq in the absence of interband 
couplingoS In our formalism, this last step amounts 
to approximating the vertices of S'lhb by local vertices 
(in time) using the atomic equations of motion for the 7| 
field. 

Since the quadratic and quartic terms of the action 
S'lhb are sufficient to recover the action of the t-J model, 
we expect higher-order vertices to vanish at the same 
level of approximation. Consider the contribution of or- 
der t/U to S'lhb due to r™j^ ^^j^ (Fig. |^a). This term 
involves the product 



J 



i^^hh (^2)7;'T (^3)7r"T (T-^)7rT (T-^TrT ir[), (3.28) 

where the fields are evaluated at different times and three 
different sites (r, r' and r"). If we now approximate this 
term by using the atomic equation of motion, we ob- 
tain the product [7r|(''')]^7r't(''')7r"T(''')[7rT(''')]^ which 
vanishes since it involves squares of Grassmann vari- 
ables. The same reasoning holds for higher-order ver- 
tices. Therefore, using the same approximations that 
lead to the t-J model, only the quadratic and quartic 
terms of the action Slhb subsist. 

We conclude that in the s trong-coupling limit, we 
can truncate the action (3.24) to quartic order in the 
fermionic fields: 



rT,r'T7r'T 



M,r'i7r'i 



S[7*, 7; - E / drj;^i9r - ^ + ^°)7rT - E / "^^^r^^^ 

r r,r' 

+ E / "^^^rl {dr-fi + U- Ai!)7rX - E / ^^7^4 
r r,r' 

"E/ '^'^(7rT^rT,r'i7r'i +C-C-) 
r,r' 

+ E / dTidT2dT3dTir\\^^;^{Ti,T2;T3, r4)7;|(Ti)7;^(T2)7ri(r4)7rt(T3). 



(3.29) 



r 



Eq. ( 3.29D is the main result of this paper. It gives 
the effective action of the Hubbard model in the strong- 
coupfing limit t/U <C 1. The fields 7-f and 7^ correspond 
to particles propagating in the LHB and the UHB. They 
are coupled by the interband hopping matrix t^i and the 
two-particle atomic vertex and interact with the 

spin fluctuations which show up in the dynamics of O. 
Note that one can describe the interaction with spin fluc- 
tuations by a SU(2) gauge field with A° its time compo- 
nent. 

The action ( |3.29 ) is suitable for a perturbative calcu- 
lation of the free energy and the LHB Green's functions. 
However, except at half-filling, it does not allow to com- 
pute the UHB Green's functions. In a hole doped system, 
closed loops of 7^ particles do not vanish (see Appendix 
^ so that the arguments used above for the LHB do not 
hold for the UHB. [Even at zeroth order in t/U, the effec- 



tive action Suhb contains an infinite number of terms.] 

An interesting aspect of this action is that the 
strong local Coulomb repulsion is already present in the 
quadratic part of the action. To u nders tand the physical 
meaning of T^^^ we rewrite Eq. ( 3.13 ) for real frequen- 
cies {iu! to) as 

EuuB — -Blhb = U + (3.30) 

where we have identified lu2 (^3) with the energy of the 
particle in the UHB (LHB). Eq. ( ^.30D suggests that we 
can interpret F^^j^ ij, the residual interaction (given that 
a part of the interaction is included in the quadratic ac- 
tion) between two particles siting at the same site. Alter- 
natively, one can view the U term in the action as a mean 
field, and T^^^ n ^ the correction to this mean field. 



Finally, we note that the action ( p. 29 ) without the 
quartic term is similar to that obtained within the 
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large-/7 Hartree-Fock approximation where the SU(2) 
spin-rotation invariance is maintained by introducing a 
fluctuating spin-quantization axis fl in the functional 
integral.Ej Omitting the quartic term (Fj^j^ -^^) is however 
in general not possible without missing processes of order 



t/U. For instance, this tert 
derivation of the t-J modeli 



,plays a central role in the 



D. Strong-coupling diagrammatic perturbation 
theory 



The action (3.29) can be used as the starting point 
for a perturbative calculation with respect to t/U. We 
discuss in the following the computation of the free en- 
ergy and the LHB Green's function to first order in t/U. 
Formally, we have 



vnz[n], 



'•S[7*,7;J^] 



(3.31) 



and 



I?ne-5["1 {RaUGl (a, b) {RI)u, (3-32) 
^T(a,&) = / mia^lhe-'^''^'''"\ (3.33) 



1 



z[n] 



where S[n] is the action for the spin degrees of freedom 
and Q^{a,b) the Green's function for a given configura- 
tion of n. Note that the Green's function of the orig- 
inal fermions c (those involved in the definition of the 
Hubbard mod el) is not merely — (7a7h) (see Appendix 

the 



Eq. ( |3.32| ) shows that the fermionic field c 

LHB can be expressed as 



,LHB 



LHB 



= (^r)aT7rT 



(3.34) 
(3.35) 



where the last line is obtained by writing the SU(2)/U(1) 
rotation matrix i?r as 



Rr 



(3.36) 



with the constraint |z^-|-| -I- |z^| | = 1. Eq. ( 3.35 ) is famil- 



iar from the Sch|«inger-boson slave-fermion formulation 
of the t-J modelO 



J 




FIG. 8. An example of self-energy diagram for the UHB 
particles which is 0(1) in t/U . 



We now introduce the propagator 

^(o)-i 



(3.37) 



To zeroth order in t/t/, the LHB Green's function is tjj 
and the partition function is given by 



(0) 



Z(o)[J^]=det(-^j°)-^). 



(3.38) 



As discussed in Sec. flH C| Q^^^ is not the UHB Green's 
function to zeroth order in t/U except at half- filling. 
Even in the absence of interband transition, there are 
0(1) corrections to ^j"'' due to the non-vanishing of 
closed loops of_7| particles in the (hole) doped system 
(see Appendix |^). A finite Oil) self-energy correction 
to Cf^ is shown in Fig. 0. 



(a) 



(b) 



(c) 



A' 



X 

>: 



-v 



FIG. 9. Diagrammatic representation of the effective ac- 
tion of the spin degrees of freedom at order 0{t/U) . (a) 5*1 [f2] . 
(b) ^([fJ]. (c) [f2]. Slashed dashed lines indicate interband 
transitions. 

To first order in t/U ^ there are three contributions to 
5[r2] shown diagrammatically in Fig. ^: 



a,b,a' ,b' 

5j[r2] = X] X] ^-^j^(a,6;a^50G|(a^a4)^tT(^4,a3)^j"^(a3,a2)^TT(^2,al)G 

a,b,a' ,b' ai,bi 

xGiib', bi)ii^{bi, b3)g\°\b3, b2)in{b2, 6i)Gx(6i, 6), 
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Snn] = - E^n,Ti(«'^;^^«0GT(6^&4)^n(^4,63)^^°^(&3,62)^Ti(^2,6l)Gi(&l,6) 

a.b.a' ,h' a-i ,bi 

xGiia', a4)f|-f (04, as)^^^ (03, 02)^11(02, ai)G^ (ai, o). 

I 



(3.39) 



The effective action S[fl] of tlic spin degrees of free- 
dom to order 0{t/U) can be computed ex actly at half- 
filling. The zeroth-order contribution [Eq. ( 3.38| )] yields 
the Berry phase term 5'^ to lowest order in~A". When 
computing S'l, S[ and 5" we may replace ^j*''' by G_u 
since the neglected terms will be of higher-order in t/C/.E3 
Furthermore, since S'l, 5^ and S" describe (virtual) in- 
terband transitions, we can ignore the effect of the Berry 



(a) 



- + --- 



--- + - 



Q 



(b) 



-.(0) 




phase term and replace G| by G| 

(irtrT,r'i^r'i,rT 



We then find 



r.r' 



Q 



2U 



r,r' 



(3.40) 



(c) 



In Eq. (3.40), the exchange interaction is assumed to 
be instantaneous, since the characteristic spin fluctua- 
tion energy J = At^ / U is much smaller than U. The 
last line of Eq. ^( ^.40 ) is obtained from irt.r'x^r'X.rt — 
(1 — rJr • r2r')/2.E3 Since S[ vanishes at half-filling, while 
S'{ turns out to be of order 0{t^ /U^), we verify that the 
effective action of the spin degrees of freedom at order 
0{t/U), Sb + Si, is nothing but the action of the (quan- 
tum) AF Heisenberg model expressed in terms of spin 
coherent states: 



FIG. 10. Diagrammatic representation of the LHB 
self-energy at order 0{t/U). (a) Ei. (b) E'l. (c) E" (the 
symmetric diagram is not shown). Slashed dashed lines indi- 
cate interband transitions. 

The self-energy S of particles in the LHB is defined 
by Q^^ = Q^'^ ^ — S, where is the propagator for a 
given configuration of $7. It has three contributions to 
first order in t/U (Fig. |l|): 



-7 • (3-41) 



Si(a,a') 
S'l (a, a') 
E"(a,a') 



b,b' 

Y ru,n(a, b; a\ b')Gi{b', h)ii^{bi,b2)g\"\b2, 53)tTi(&3, b4)Gi{bA, 5), 

bM',bi 

^ ^ iTi(«, ai)Gi(ai, a2)Fi/^^^^(6, as; a^ 50Gx(6^ 5i)£n(&i, ^'2)af ^(62, 63)t^T(&3, 64)Gt(64, 6) 

b,b' ai,bi 

+ EE^UJi(«'^;^^«2)GT(6^5l)^TT(6l,62)^j"^(62,63)^u(^3,64)Gx(64,fo)Gx(a2,al)^iT(al,a'). (3.42) 

b^b' aijbi 



Si describes indirect hopping processes between next- 
nearest-neighbor sites which occur via virtual transitions 
to the UHB. In the t-J model, these processes show up 
in the pair-hopping term. 

So far our discussion has been confined to a formal 
framework. For practical calculations, further approxi- 
mations will be necessary. Whereas the fermionic fields 
7o- can be integrated out in a systematic t/U expansion, 
one still has to carry out the functional integral over the 



unit vector field O. 

The simplest approach consists in expanding around 
a broken-symmetry ground-state by making a saddle 
point approximation on the spin variables Hr- Differ- 
ent choices, corresponding to AF, ferromagnetic or spiral 
orders, are possible. This approach, which is discussed at 
length in Ref. can be justified by taking a "large-5"' 
semiclassical limit of the Hubbard model. The spin-^ co- 
herent states can be promoted to spin-S coherent states 
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(with S arbitrary) by w riting the SU(2)/U(1) rotation 
matrix Rr- as in Eq. (Uej) and generahzing the constraint 
l-^rtP + kriP = 1 to l^rTp + Im P = 25*. The latter 
equahty allows to write Zrt = V2^cos(6'r/2)e"4('^'-+'^'-' 
and = V^S sm{dr/2)ei^'^''~^'''> where the choice of V'r 
is free and corresponds to the U(l) gauge freedom. The 
action of the Hubbard model is obtained from ( pS| ) with 
the replacement A° 2SA° and i 2Si. [Note that 
the creation of a (fermionic) hole corresponds to a total 
removal of the local moment.] In the limit S —^ oo (with 
tS = constant), the Berry phase term suppresses quan- 
tum fluctuations of CI. The spin variables become classi- 
cal and do not fluctuate at zero temperature. The leading 
correction in 1/5 gives the spin- wave modes around the 
broken-symmetry ground-state. 

At half-filling, one can also try to c ompu te G by di- 
rectly averaging Q with 5[r2] [see Eq. (3.32)], since the 
effective action of the spin degrees otireedom, S'Hcis[J^], 
can be analyzed by various methods. l3 Since S'Hcis[J^] is 
not a Gaussian action, this seems a priori a very difficult 
task. We show in Ref. ^ how this technical difficulty 
can be circumvented by using the Schwinger-boson mean- 
field theory of the Heiscnberg model. The advantage of 
this approach is that one starts from a good description 
of the magnetic properties of the system (which are de- 
scribed by S'Heis): absence of long-range order at finite 
temperature, exponential divergence in 1/T of the mag- 
netic correlation length, etc. Moreover, this formulation 
is a good starting jppint to study how spin fluctuations 
affect hole motion.E3 

It should be pointed out that Eqs. ( 3.32| ) and ( 3.33| ) 
cast the Hubbard model in the form of a spin-fermion 
model, which can be seen as the strong-coupling coun- 
terpart of the "weak-couoling" spin-fermion model that 
has been studied so far.t2l There are however important 
differences between the latter and our formalism: (i) In 
the "weak-coupling" spin-fermion model, the effective ac- 
tion of the spin degrees of freedom is generally assumed 
to be Gaussian. In our approach, the effective action 
S[fl] is obtained exactly (to order 0{t/U)) by integrat- 
ing out the fermionic degrees of freedom, (ii) The spin- 
fermion model omits coupling between spin fluctuations 
and charge degrees of freedom. Fermions interact with 
spin fluctuations via a spin-spin interaction. In our ap- 
proach, there is an intimate coupling between charge de- 
grees of freedom (i.e. hole motion) and spin fluctuations, 
which shows up in the f2-dependence of the intersite hop- 
ping matrix, (iii) The spin-fermion model is restricted to 
the weak-coupling regime, whereas our -approach is ap- 
propriate to the strong-coupling regime.EB 

Finally, we would like to discuss the relatio n to t he t- 
J model. Starting from the effective action (3.29) and 
integrating out the UHB, one recovers the action of ths. 
t-J model in the spin-hole coherent-state path integral.E3 
A further change of variables yields tlie, action of the t-J 
model in the slave fermion formalism.EZl [See Sec. I II.B.3 
of Ref. |9[] The connection between the action ( 3.29| ) 



and the t-J model is however not straightforward and 
involves some subtleties regarding the proper time or- 
dering in the functional integral.Ej As a result, both for- 
mulations should be considered as two different descrip- 
tions of the strong-coupling limit of the Hubbard model. 
For instance, indirect hopping processes between next- 
ncarc st ne ighbors appear in the Gaussian part of the ac- 
tion (3.29), while they are described by the pair- hopping 
term in the t-J model (this term involves a three-body 
interaction in the standard formulation of the t-J model). 
But the most interesting aspect of our formulation is the 
use of spin-particle-hole coherent states. By clearly dis- 
tinguishing between charge and spin degrees of freedom, 
spin-particle-hole coherent states open up new possibili- 
ties for analvgjng hole motion in presence of strong spin 
fluctuations .c3 



IV. SUMMARY AND CONCLUSION 

There are two main ingredients in the derivation of the 
Heisenberg model from the Hubbard model at half-filling: 
First, one should deal with spin operators (instead of the 
original fermionic operators) in order to account for the 
local moments that appear in the Mott insulating phase. 
Then, the effective Hamiltonian for these spiftjOperators 
is derived by expanding with respect to t/U£h 

In this paper, we have proposed a generalization of 
this procedure to the doped case. Our approach is based 
on the introduction of spin-particle-hole coherent states 
which generalize the spin-i coherent states by allowing 
the creation of a hole or an additional particle. They also 
generalize the spin-hole coherent states introduced in the 
context of the t-J model.o The spin-particle-hole coher- 
ent states can be used to derive a path integral and natu- 
rally lead to a spin-rotation-invariant slave-boson formu- 
lation of the Hubbard model (Sec. ||). By performing two 
successive Hubbard-Stratonovich transformations of the 
intersite hopping term, the path integral can be recast in 
a from which is suitable for a perturbation expansion in 
t/U . The action S'[7*, 7; ^2] can be expressed in terms of 
two Grassmann variables (7f ,7x) and a unit vector field. 
A singly occupied site is described by a spin-i coherent 
state irij.), while 7^ (7^) describes a particle propagating 
in the LHB (UHB). The fermionic action has interac- 
tion t erms to all orders, given by the atomic vertices F^ 
[Eq. (3.24)]. The fermions also interact with spin fluc- 
tuations which show up in the dynamics of fi. In the 
strong-coupling limit t <^U , the action can be truncated 
to quartic order in the fermionic field [Eq. ( 3.29| )] and 
used as the starting point of a diagrammatic expansion. 
At half-filling and to order 0{t/U), we recover the action 
of the (quantum) Heisenberg model. 

Since the fermionic field can be systematically inte- 
grated out within a t/U expansion, the main practical 
difficulty comes from the dynamics of il (spin fluctu- 
ations). The simplest approach consists in expanding 
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about a broken-symmetry ground-state my making a sad- 
dle point approximation on the spin variables Sir- Differ- 
ent choices, corresponding to AF, ferromagnetic, or spiral 
orders are possible. This approach is discussed in detail 
in Ref. The spin-particle-hole coherent-state path in- 
tegral turns out to be a very convenient tool for studying 
spin waves about a broken-symmetry ground-state in the 
strong-coupling limit. 

The main characteristic of the spin-particle-hole 



coherent-state path integral lies in the clear distinction 
between charge (7) and spin ($1) degrees of freedom. This 
opens up new possibilities for the computation of the 
Green's functions. One can first calculate the Green's 
functions for a given (time-dependent) configuration of 
H and then perform the average with the effect ive ac tion 
S[fl] of the spi n degrees of freedom [see Sec. HID and 
Eqs. (3.3S-3.33)]. This program is carried out in Ref. E|. 



J 



APPENDIX A: PATH INTEGRAL FOR THE SINGLE-SITE HUBBARD MODEL 



Introducing (M — 1) times the resolution of the identity (2.13) in Eq. (2.22), one obtains 

\fe=l ^ / 



(Al) 



with e = 13/M and the boundary conditions: fJo = J^Af, em = sq, pu = Po, dM = do, and f^M = —fao- The 
projection operator Pk is defined by Eq. (2.15) with the spin-quantization axis determined by the imit vector Jl^. It 
can be written as 



Pk = 



k _ — ieA 



2tt 



n 



27r 



(A2) 



Since the constraints Q^^-* and Q^^"* commute with the Hamiltonian H — /iiV, we obtain 



in the limit e — > 0, where 



fc=i / fe=i 

M ,„ \ A/ 



Kk=l 



(A3) 
(A4) 



dXk = 



k=l 



2tt 



n 



€dX\ 



2Tr ' 



kk = H-nN + iX^^'^Qi'^ + I xf^Q 

(T 

{nkXk\Kk\^k-iXk-i) 



(2) 
a 5 



k.k-1 — 



(A5) 



Using the expression of the scalar product [Eq. ( 2.12 )], the partition function is finally written as 

^ / A/ ,„ \ r Af 

^at = AA^' / n -^dQdCkdXk expi ^ [-a\Ck? + CkCk-i + /ffcPl./Tfe-iW-i((f^fc|f^fe-i) - 1) - 

\k=l ^ ) U=l 



(A6) 



Taking ae = — 1 and neglecting the overall normalization constant M^^ (see the discussion in section |ll B| ), 
there is now no difficulty to take the continuum time limit. Kk.k-i — Kk^k is evaluated by writing the Hamiltonian 
as 



H-t,N = Ud^d~fiJ2f^J-^ 



(A7) 
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where the spin-quantization axis flk is chosen. The constraints Q*-^-* and Qi^"* being also defined with respect to the 
spin-quantization axis fl^., we obtain 

Kk,k = U\dk\^ - A*E + ^^k\\^k\' + \Pk\' + \dk\' - 1) 

+*AS(/ffc/T/c - \Pk\' - |4n + ^xf,!{flJlk - |4n. (A8) 
We also note that (r2fe|r2fc_i) — 1 = (r2fe|r2fc_i) — (Ofe|r2fe) yields a term — (OjO) in the continuum time limi t. He re 



the dot denotes a time derivative. The final expression of the action ^at in the atomic limit is given by Eqs. (2.2E) 



APPENDIX B: EVALUATION OF {(IXICrCr'I^X) 
In this appendix, we prove Eq. (|2.3l| ) (we do not write explicitly the time index k). We introduce the operator 

$1 = 5](i?r).'.cL, (Bl) 
a' 

which creates a particle with spin a in the spin reference frame determined by ft^. In the space spanned by B'^^ |, 
we have 4>ra — %a-- We then deduce 

{n, C\clcr' \n, = ^(f^, C\jl{RlRr')aa'%'a'\^, C) 

= jlRlRr'Jr'{nX\n,0, (B2) 

where 7r| = e*Pr/rT ^^^d 7rx — p*c?r/rX- The intersite term in the action ( |2.34 ) follows from Eq. (B2). 



APPENDIX C: GREEN'S FUNCTIONS IN THE ATOMIC LIMIT 

In this section we compute the connected atomic Green's functions 

Gfalb,} = i-lf'iTrla, ' --la^ll ' ■■ll).t,c, (CI) 

written here in the operator formalism (T^- is the imaginary-time ordering operator). We consider a single site and 
drop the site index. 

Let us first consider the Green's functions determined by . In the basis B'^, the operators and the Hamiltonian 
H — jiN can be written as (in matrix form) 

7T = 

/ 

H - fiN = \ -n I . (C2) 
\0 U-2n 

We easily find the time dependence of the operators ^a\T) = U i—T)j!^'' U (t) (where U{t) ~ e^T{H-nN) -j-j^g 
evolution operator): 

7T(r) = e^^7T, l\{r) = ^\ , 

7i(T) = e(^-^)^7i, 7l(r) = e-(^-^)^7l- (C3) 
The single-particle Green's function Gct"''(t) = — (7V7(r)7j5.(r)) is found to be 






1 
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7(0) 



Z. 



(0) 



re 



(C4) 



The final expressions are obtained in the low-temperature limit and we have used Z^^^ = 1 + e^^^ + e'''^^^ ~ e'^'^. 
The latter approximation neglects the contribution of the empty and doubly occupied sites, in agreement with the 
discussion of Sec. II B. As expected, zf^ is the partition function for a given spin direction. In Fourier space, we 
obtain 



Gf\iuj) = [iu + ii-U)-'. 
In the same way, we can calculate the two-particle Green's function. For instance, we have 



(C5) 



X [0(t4 - T2)e{T2 - T3)0{t3 - n) - 0{t3 - T2)0{t2 ~ Ti)0{n - Tl) 



(C6) 



where Tr[{/(/3)(7|7i-)^] ~ We do not distinguish between G^^*-*^-* and G^\ since the effect of the Berry phase term 
on G^'^ (> 2) can be ignored (see Sec. [II). One can verify that G^^ ||(ti, r2; T3, T4) coincides with its disconnected 
part 



TT,TT(^i'^2;r3,r4) - G'"'{t, ~ T,)Gf\T2 - T4) - Gf'(ri - T,)Gf\T2 - ra) 



'[^(-Tl + T^)0{-T2 + Ti) - 0(-ri + Ti)0{-T2 + T3)], 

so that Gji^ii = 0. Similarly G^^^ = 0. By using the definition of the third cumulant 



(C7) 



aia2a3,bib2b3 



= G 



III 

010203, f)i62')3 



0161^0203,6263 



Mo) ^iic 
'-^0262 0103,6163 



(0) ^IIc 



- Gi"i G 



0363 '^102, 616; 



(0) ^IIc 



+ G 



0261 ^0103,6263 



_^(0) ^IIc I ^(0) ^IIc I ^(0) ^IIc I ^(0) ^IIc _ p(0) pile 

0361 0102,6263 ' 0362 0102,6163 0162 0203,6163 ' 0263 0103,6162 0163 0203,6162 

_r'-0) r^(o) ^(0) _j_^(o) ^(o) ^(0) , ^(o) ^(o) ^(0) 

'-^Oi6i^0262'-^0363 "t" "-^0162 '-'0261 "^0363 ^0i6i'-'0263'-^0362 



MO) MO) r^(o) 

'^Oi62^0263'-^036i 



MO) MO) MO) ,Mo) MO) Mo) 

"-"0163*^0261 "-^0362 '-^0i63'-^0262"-'036i' 



(C8) 



one also obtains GiH^L 



'aaa.crcrcr — 0- [Note that this verification requires to consider 6!=360 different time sectors.] It is 
clear that this result holds to all orders, i.e. 



(C9) 



The Green's function G^^-^^ can be calculated in the same way: 

Gj\,^^(ri,r2;T3,r4) = 6^.(^1 -3)+(M-C/)(r2-r4) _ ^^^^(^^ _ ^^)^(^^ _ 

-'0{r2 - n)e{T3 - Tl) [0in - T2) + 0{t^ - T3)] + e'3(M-C/)0(^4 _ ^3)0(^3 _ r^)e(^r^ _r2)]. (CIO) 

In Fourier space, we obtain from Eq. ( |C10 ) the connected part 

G^^^_^^{iuji,iuj2;iuj'i{i^4)) = -G^^\iuJi)G^^\iuj2)G^^\iujz)G''^\iuj4){iuj2 - ii^s - U) 

-^^^Gf\iu,)Gf\iu:,) ~ eP^^-^)Gf{iu2)Gf{iu^A)\ , (Cll) 



1 



ibJi + iui2 + 2fi — U 
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where UJ4 = uji + uj2 — 1^3 is fixed by energy conservation. G^^'^ is related to the two-particle vertex by 



|l(iwi,icj2;«t^3(«^4)) = -Gf\iuji)G'--l^\iu!2)Tii^ii{iuji,iu!2;iuj3, {iuj4))Gf\iuj3)Gf\iuj4) 



.(0) 



From Eq. (Cll), we then deduce 



(C12) 



(C13) 



in the low-temperature limit, i.e by neglecting the second and third terms of the rhs of Eq. (Cll). These t erms 
correspond to the contributions of the empty and doubly occupied sites, which, according to the discussion of Sec. II B , 
should be discarded. 

Let us now consider the effect of the Berry phase term 



(C14) 



on the single particle Green's function G. The latter can be calculated from 

1 



Gi(r) = - 



7(0) 

1 



dX / I?[/,e,p,d]e*(T)p(T)/T(T)/f(OK(0)e(0)e 



Z. 



(0) 



dX / I?[/,e,p,d]d(T)p*(r)/i(T)/nOMO)d*(0)e- 



(C15) 



by integrating explicitly the bosonic (e,p, d) and fermionic (/o-) fields and the Lagrange multipliers. The calculation 
can be found in Appendix C of Ref. p9. Neglecting the Berry phase term, one recovers the results obtained above 
[Eqs. (C4)]. The first order correction in A'^ = (r2|ii) is 



g\'\t) = J dnG^{'\T-n)A"in)Gf\n), 
g\'\t)=- f dnGf\r-rM°iTi)Gf\r,) 



(C16) 



in the limit T ^ 0. These results are easily extrapolated to higher orders as (in matrix form) 

G^ = Gf^ + Gf^A'Gf + GfA^Gf^A^Gf^ + ■■■ 



Gx = Gf - G^A^Gf^ + Gf^A^G^A'^Gf^ 



or, equivalently, 



g; 



Gf^-'~A^, 



GJ^ = Gf^-' + A'>. 

Note that these Green's functions, like the gauge field A'^, depend on the site which is considered. 



(C17) 



(C18) 



APPENDIX D: CLOSED LOOPS IN THE STRONG-COUPLING PERTURBATIVE EXPANSION 

In this Appendix, we discuss the vanishing of closed loops of 7^ particles in the strong-coupling perturbative 
expansion. 

Since the Green's functions G^'^ are related to the vertices F^, a diagram can always be expressed in terms of the 
G^'^'s. Using Eqs. (|C3|), one can show that G^'^ is a sum of terms proportional to (—1)''^ Y[f=i {ti — Tp(i)) where P 
is a permutation of [1, B\ satisfying ap(i) = (i S [1, i?]). Each closed loop of 70. particles will then give terms of the 
type Gct(to — Ti)taa ■ ■ ■ taaGa{Tn — tq) where a is the (conserved) spin along the loop and the integer n varies from 2 
to 00. Terms with 71 = 1 correspond to anomalous diagrams and should be discarded. By summing over the different 
values of n, we obtain Q'^\r,v]T = 0^) — Gcr(r,r;r = 0^) where q'^'' is defined by Eq. (3.37). In the absence of 
particles in the UHB, both ^|''-'(r,r;T = 0^) and G|(r,r;T = 0~) are equal to zero. We conclude that closed loops of 

7l particles vanish. On the other hand, since Q^^^ (r, r; r = 0^) — G| (r, r; r = 0^) = n—1, closed loops of 7^ particles 
do not vanish, except when the density of particles n equals 1 (half- filled case). 
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APPENDIX E: GREEN'S FUNCTION OF THE "PHYSICAL" FERMIONS 



In this section, we relate the "physical" Green's function — (cacj) to the propagator of the 7 field. Consider the 



fermionic field (j) defined in the spin reference frame determined by O [Eq. (Bl)]: 



r J (7(7' H^va' 



(El) 



Since = 7 (in the space spanned by B'^^ ) , we deduce 



(cr.(r)c:,,(r')) ^-Y. ((^r(r))..,7r..(r)7r%.(r')(4(^')). 



(E2) 



For fermions in the LHB, we then find 



^LHB 



(rr,r'r') = -((i?r(r)).T7rT(r)7r* t(^')(4 (^'))t.>- 



(E3) 



Here 7^ = e*p/f , 7^ = and the mean values are to be taken with the action S[f, e,p, d, A; O] [Eq. ( ^.34| )]. The 

Green's function C/^hb ^^^^ expressed as the functional derivative of the partition function calculated in the pre sence 
of ex terna l sources. It is then possible to carry out the steps leading to the effective action S'[7*,7;r2] [Eqs. ( 3.21 ) 
and ( ^.24 )]. Taking the functional derivative, one finds that Q^^^ is given by Eq. ( pil3[ ) where the 7 field is now an 
elementary field and the mean value has to be taken with the action £'[7*, 7; n]. 
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